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Abstract. We use G-band and Ca n H image sequences from the Dutch Open Telescope (DOT) to study magnetic elements 
that appear as bright points in internetwork parts of the quiet solar photosphere and chromosphere. We find that many of these 
bright points appear recurrently with varying intensity and horizontal motion within longer-lived magnetic patches. We develop 
an algorithm for detection of the patches and find that all patches identified last much longer than the granulation. The patches 
outline cell patterns on mesogranular scales, indicating that magnetic flux tubes are advected by granular flows to mesogranular 
boundaries. Statistical analysis of the emergence and disappearance of the patches points to an average patch lifetime as long 
as 530 ± 50 min (about nine hours), which suggests that the magnetic elements constituting strong internetwork fields are not 
generated by a local turbulent dynamo. 
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1. Introduction 

In this paper, we address the appearance and lifetime of mag- 
netic elements that intermittently show up as bright points in 
the internetwork areas of the quiet sun. The context is the nature 
of quiet-sun magnetism, its dynamical coupling to transition- 
region and coronal fields, and the existence of a local turbulent 
dynamo. 

The more familiar network bright points are similar They 
constitute the magnetic network which partially outlines the 
boundaries of supergranular cells, and have long been recog- 
nized to represent strong-field magnetic elements, which are 
traditionally mode led as flux tubes. They we re first observed 
as magnetic kno ts ( Beckers & Schrote3ll968 ) and as "filigree" 
jPunn & Zirkei^n ^73l) r esolve d into strings of adjacent bright 
points by Mehltr etteil (Il974l) . Prompted by Muller (1977), 
Wilsonl (119811) showed that faculae, filigree, and bright points 
in wide-band Ca n H filterg rams are manifestations of the 
same phenomenon. MulleiTl 1983|) subsequently introduced 
the name "network bright point" (NBP) and initiated an 
extensive hterature obser ving them as G-band bright points 
dMuUer & Roudieij Il984l) . In particular the G-band studies 
with the f ormer and the pres e nt Swe di sh solar telescope on 
La Palma (iB erger et al."l995', 'l998blR '2004 iBerger & Title 
19961 laOOlFlLofdahl et al. 1998; van Ballegooijen et al.lll998t 
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Wiehr et al.l 120041 iRouppe van der Voort et al.l l2005h estab- 
lished that NBPs are brightness manifestations of the small, 
strong-fiel d magnetic elements that make up the magnetic 
netwo rk ( Chapman & SheelevI 19681: Livingston & Harvey 



1969^ 



11972 



Howard & Stenflo 
Stenflol Il973h . These have been 
flux tubes since the 



Frazier & Stenflol 11972 : 
modeled as magnet ostatic 
pio n eering work of .Spruij ( 1976 , 
1977h inspired by IZwaanI (Il967h . Their hot- wall e xplana 



tion o f the photospheric b rightness enhancement dSpruit 
19761; ISpruit & ZwaaiJ 11981 ) wa s recently v erified by the 



MHD simulations of iKefler et al 
whic h crown a long 



(120041) . 

ing (e.g., iKnoUcer & Schiisslei 



(120041) and ICarlsson et al 
effort in flux tube model- 



Solanki & Briglievic 



1992t 



19881; Keller et al 



Grossmann-Doerth et al.l 



199C; 



1994 



19981; ISteiner et al.lll998l:ISteineijl2005l) . NBPs serve as tracers 
of strong-field flux tubes, especiaUy in the Fraunhofer G band 
(CH lines around 430.5 nm) and the CN band at 388 nm (cf. 



Rutten et al.ll2001h . but with the caveat that not all magnetic 



featur es produce observable bright points (Berger 
200lh . 



The nature of solar magnetism in the quiet-sun inter- 
network areas, i.e., the supergranular cell interiors bordered 
incompletely by filigree chains and NBPs that form the 
network, is less well est abhshed, but is presently under 



intense scrutiny (e.g . . iDominguez Cerdefia et al 



Sanchez Almeida et al 



Lites & Socas-Navarro 



2003t 
2004t 



Socas-Navarro et al 



Socas-Navarro & Lites 



2003 



2004 



2004 
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Fig. 1. A sample image pair from the Dutch Open Telescope (DOT). Left: G-band image. Right: co-temporal and co-spatial Ca ii H 
image, clipped to improve contrast. The network is readily identified in the Can H image as regions with enhanced brightness. 
In the G-band image these areas contain many tiny bright points, often arranged in strings (filigree) within intergranular lanes. 
These are NBPs. Close inspection shows the presence of isolated IBPs in the internetwork regions. 



Truiillo Bueno et al.l 120041: iManso Sainz et alj |2004 . A full 
spectrum of field strengths seems to be ubiquitously present 
in the internetwork at small spatial scales, with the stronger 
elements residing in intergranular lanes. When strong enough, 
such internetwork elements are sufficiently evacuated in the 
low photosphere to appear as internetwork bright points (IBPs) 
that are quite similar to NB Ps, but more i solated. Their exis- 
tence was already noted by Mullei ( 1983 b. but only recently 
a more detailed study of photos pheric IBPs was presented 
by ISanchez Almeida et alj (2004). They report a density of 
0.3 IBP per Mm^ and lifetimes of a few minutes, shorter than 



the aver age NBP lifetime of 9.3 min measured bv lBerger et al 
(Il998bl) . 



It is advantageous to combine photospheric imaging in the 
G band with simultaneous co-spatial imaging in the core of the 
Can H line to study IBP appearance, patterning, and lifetimes. 
In Ca 11 H & K images sampling the low chromosphere, both 
NBPs and IBPs show up with a larger brightness enhancement 
over the surrounding area than they do in G-band or contin 



uum images sampling the photosphere (cf. Fig. 2 of lLites et al 



1999 ). Ca n H imaging therefore provides a better diagnostic to 
detect isolated IBPs, for which contrast is more important than 
sharpness. The latter is lower for Can H&K images due to 
strong scattering, and possibly also due to increasing flux tube 
width with height. G-band IBPs vanish sooner than their Ca ii 
counterparts when they diminish in brightness. G-band imag- 
ing at the highest resolution is needed to faithfully render the 
intricate brightness structure of individual magnetic elements 
dBerger et alj |2004|) . but synchronous Can imaging provides 



better location detection especially for isolated ones, as IBPs 
often are. 



Intermittent IBPs o bserved in C an H&K 
"persistent flashers" by iBrandt et al 



were called 

d 19921 1 1994 . who de- 



scribed a particular Can K2V grain which appeared and dis- 
appeared during multiple hours while migrating from the cen- 
ter to the boundary of a cell. The grain followed a flow 
path determined from independent granulation tracking as if 
it were a cork floating on the solar surface. These flashers 
figured in the extensive debate whether all briefly-appearing 
bright grains in Can H&K image sequences, in particular 
at the H2V and Kay off-center wavelengths, represent non- 
magnetic acoustic shocks o r mag netisni-constrained phenom- 



ena (e .g., Sivaraman & L ivingston 1982': 'Rutten & Uitenbr oekl 
I1991I; I von Ue xkuell&Kneer. .1995i: jCarlsson & Stein, , 1997 
Nindos & Zirin| 
(Sivaraman et al. 



19981: IWorden et all [1999; Lites et al] Il999 ; 



2000). The conclusion is that the acoustic in- 



ternetwork grains appear only a few times with a modula- 
tion of roughly three minutes, whereas the strongest internet- 
work magnetic elements stand out by their "persistent flasher" 
character: they possess an apparent location memory which 
significantly exceeds a few three-minute cycles. Many such 
longer-lived (while modulated) internetwork brightness fea- 
tur es were noted in ultr aviolet image sequences from TRACE 
by iKrijger et alJ (120011) . who displayed one in their Fig. 7. 
A higher-resolution exampl e from the Dutch Op en Telescope 
(DOT) is shown in Fig. 2 of iRutten et al] (l2004al) . 



In this paper, we study IBPs in synchronous G-band and 
Ca n H image sequences with high resolution and fast cadence. 
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Fig. 2. Example of a magnetic patch, visualized by partial Ca ii H images (x-y cutouts, top row), the associated Ca ii H x-t slices 
(middle row), and the corresponding x-t slices from the IBP map sequence (bottom row). The images shown here are a selection 
of the data at intervals of At = 340 s starting at f = 240 s and Ay = 0.426" starting aty- 0.568", respectively. The rightmost 
panels are averages over all data collected in the duration of the sequence (top panel) or in the y interval shown here (lower 
panels). The y location of each x-t slice is shown by a dashed line in the associated x-y panel, while the time of the latter is 
shown by a dashed line in the corresponding x-t shce. 



as a sequel to the remark by Rutten et al. ( 20043) that the IBP 
identified as magnetic flasher there persisted over the full se- 
quence duration of 54 min. The main issue is whether briefly 
appearing IPBs systematically portray longer- lived flux con- 
centrations that vary temporally in their morphology and inten- 
sity. Our Ca II H imaging and statistical analysis below indeed 
suggest that this is the case. This result is important in the con- 
text of fi eld gene ration by turbulent dynamos (Catt aneolll999t 
lEmonet & Cattaneol 120011; ICattaneo etalj|2003h . of the field 
topology surrounding network at chromospheric heights (cf. 
[Schrijver & Title 2003), and of the internetwork contribution to 
coronal heating by wave dis sipation (e.g., iHevvaerts & Priesli 
1983h and reconnection (e.g., |Parkeilll988h imposed by photo- 



spheric foot point motions and topology evolution. 

2. Observations, data reduction, and patch 
identification 

We use a double image sequence of a quiet area at the disk 
center recorded by the DOT from 8:40 to 9:39 UT on June 16, 
2003. The sequences consist of speckle-reconstructed images 
taken at a 20-s cadence in the Fraunhofer G band with a 
10- A filter centered at 4305 A, and synchronous, speckle- 
reconstructed images taken in the Ca ii H line (3968 A) with a 
narrow -band filter (FWHM 1.3 A) at line center A sample pair 
of G-band and Can H images is shown in Fig. [T] Details on 



the telescope, its tomographic multi-wavelength imaging and 
image acquisition, and the speckl e reconstruction and standard 
reduction procedures are given in lRutten et al. ( 2004bl) . 

The images were carefully aligned and destretched us- 
ing Fourier correlation. After clipping to the common field of 
view (68 X 63 arcsec^), the resulting sequences consist of 178 
speckle-reconstructed image pairs of 962 x 894 square pixels of 
0.071" size. Each sequence was cone-filtered in Fourier space 
to remove features that travel with apparent horizontal speed 
exceeding the 7 km s ' sound speed. The image sequences may 
be downloaded as movies from the DOT databas^l together 
with synchronous blue and red continuum sequences not used 
in this analysis. 

In order to visually search for recurrent IBPs, we employed 
an interactive three-dimensional "cube sheer", which dissects 
both data cubes simultaneously in x-y, x-t, and y-t slices with 
continuous (x, y, t) selection controlled by mouse movement. 
For example, IBPs show up in the x-t slice only when not drift- 
ing in y, but a slight wiggling of the y coordinate then helps to 
track the three-dimensional "world-line" of the feature through 
both data cubes. We so found many IBPs that are intermittently 
present while migrating slightly in x and/or y. They are easiest 
to detect in the internetwork parts of the Ca ii H sequence, but 
are accompanied, at least for part of the time, by far sharper co- 
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spatial G-band IBPs in the underlying photosphere. Often, we 
find groups of recurrent IBPs that appear associated through a 
shared magnetic structure providing a longer-term spatial lo- 
cation memory. We call such linked IBP groups "magnetic 
patches". They frequently consist of multiple strings of IBPs 
that split and merge. 

Figure |2] displays an illustrative example. The panel lay- 
out mimics cube slicing in the form of a sequence of small 
successive x-y image cutouts (top row) and a sequence of x-t 
slice cutouts stepping progressively in their y sampling loca- 
tion (middle row). The rightmost column of panels are averages 
over the plotted t and y ranges. We discuss two representative 
cases, (i) The first one is the bright IBP near y x 29.5 arcsec 
in the third x-y panel (f 15 min), which is sampled by the 
corresponding x-t slice. The latter indicates an IBP lifetime 
of about 7 min, during which it migrated leftward at about 
0.15 arcsec min"', or 1.8 kms"'. It seems to be fairly isolated, 
except that some intermittent brightness appears later at x = 5- 
6 arcsec in the same x-t slice. The next three x-t slices show a 
bright structure near f « 40 min that migrates toward larger x. 
(ii) The second example is the bright grain near y ^ 31 arcsec 
in the 7th x-y panel for f w 38 min. Its x-t slice indicates 
that it appeared abruptly at this time and then Uved for fifteen 
minutes, but the surrounding slices show enhanced brightness 
nearby in y before and after as well. True cube slicing confirms 
that the first example is indeed a continuously present magnetic 
structure which merges at f w 45 min with the second example. 

The rightmost panels of Fig. |2]encompass the third dimen- 
sion through integration over the full sequence duration (top 
panel) or over the y range shown (middle panel). The nearly 
continuous brightness in the interval x = 5-7 arcsec suggests 
that many IBPs in this area belong to a common magnetic patch 
throughout the Can H sequence duration. Fig. |4] confirms that 
this is indeed the case with the friends-of-friends patch defini- 
tion described below. Note that the IBPs of cases (i) and (ii) 
apparently intersect in this panel at t x 12 min, whereas they 
are in fact disjunct in y (as is shown in detail by the x-y panels 
in Fig. [2] and by the x-integrated y-f slice in Fig.|4]i. 

We have developed a detection algorithm to locate mag- 
netic patches made up by IBPs in order to quantify our visual 
impressions from cube slicing. It includes handling of splits 
and mergers, such as the ones discussed above. First, a mask 
to block off the network was constructed from the Ca ii H se- 
quence by temporal averaging of the entire one-hour Ca n H 
sequence, followed by 50-pixel boxcar smoothing and thresh- 
olding at the mean value, passing only lower values. The re- 
sulting mask is shown in Fig. |3] The selection is conservative 
in that medium-bright ar eas around bright n etw ork (the "in- 
termed iate" pixel class of iKrijger et al.ll200ll and lRutten et al. 
2004al) are also rejected. 

We next employed a multi-step procedure to locate IBPs 
in the masked-ofF sequences. Each image was first convolved 
with a suitable kernel to increase the contrast of small round 
features. The kernel has the form 
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Fig. 3. The temporal average of the Ca ii H image sequence 
over its one-hour duration with the internetwork mask contours 
overlaid in white. The Can H image intensity was scaled log- 
arithmically in order to show contrast in the internetwork. The 
internetwork mask is computed by taking a 50-pixel boxcar av- 
erage and thresholding at the mean value. 



where k is a chosen such that the spatial average {Knir)) - 0. 
For the G-band sequence we used a kernel with a 5-pixel diam- 
eter {D - 4). We chose a 7-pixel kernel {D - 6) for the Ca ii H 
images, in which IBPs are larger. We subsequently produced 
binary maps of IBP-candidates by thresholding the convolved 
images at a suitable level. This threshold level must be low, 
because IBPs occasionally become very weak between periods 
of enhanced brightness. The resulting maps therefore are quite 
noisy and require further processing. 

The binary map sequences were improved further through 
spatial and temporal erosion-dilation processing. A spatial ero- 
sion operation tests the local nature of a candidate IBP in a 
binary map by discarding those pixels whose surroundings do 
not match a given kernel. A dilation operation does the oppo- 
site by adding pixels around the candidate IBP pixels. For the 
Can H maps, we chose a 3 x 3-pixel kernel for the erosion as 
well as the dilation operation, so that only candidate IBPs of at 
least this size pass the test. This spatial processing was omitted 
for the G-band maps, because G-band IBPs are often smaller 
than this kernel. However, temporal erosion-dilation process- 
ing was done on both binary-map sequences using a kernel of 
three time steps (60 s) to remove short-lived features. 

Although the erosion-dilation processing removes much 
noise, there remain structures in the bina ry IBP maps that 
are not IBPs. Reversed granulation (cf. Rutten et al.ll2004a) in 
particular produces arc-shaped structures in the binary Can H 
maps that the above processing fails to remove. Many of these 
are short-lived. We therefore discarded all structures with life- 
times less than 80 s (4 images). 
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Fig. 4. Integrated x-y cutout, and x-t and y-t slices correspond- 
ing to those in Fig.|2]of the binary IBP map sequences for the 
Can H (left) and G band (right), showing only those IBPs that 
group into the central patch. Careful comparison with the right- 
most panel in the bottom row of Fig. |2]shows that several IBPs 
present there are part of another patch, such as at the leftmost 
edge of the cutout around t ^ 25 min. Each panel displays the 
integrated IBP map sequence over the third coordinate, i.e., the 
x-y map sequence cutouts integrate over time, the x-t slices 
over y, and the y-t slices over x. 



In order to retain only features with a small spatial extent, 
we keep only those features whose average maximum instan- 
taneous size in X or y expressed in units of 0.071" pixels is 
smaller than the feature lifetime measured in units of 20-s sam- 
pling intervals. Finally, we inspected all candidate IBPs visu- 
ally, either accepting or discarding them. We believe that this 
last, laborious step removes most, if not all of the misidentifi- 
cations from our sample. 

The third row in Fig. |2] shows the corresponding x-t slices 
through the resulting binary Ca n H IBP map sequence. Much 
of the brightness pattern in the rightmost slice in the middle row 
survives the IBP selection. The IBPs of case (i) and (ii) are, of 
course, properly detected by our algorithm. The apparent inter- 
section at f 12 min also remains. However, some IBPs iden- 
tified by our algorithm, for example those around x 5 arcsec 
and f w 20 min, are not visible in the y-averaged Ca ii H bright- 
ness. 

The next step is to group the IBPs into magnetic patches. 
To this end, we apply a friends-of-friends algorithm. Two IBPs 
are considered friends if their minimum separation is less than 
0.71" (10 pixels), with disregard of temporal seperation. A 
patch then consists of a group of befriended IBPs that have 
no friends outside the group. It may contain IBPs that are not 
direct friends if a string of IBPs form a path connecting them. 
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Fig. 5. The average Ca ii H image with IBPs identified in the 
Ca II H sequence overlaid in white. The IBPs appear to group in 
patches that outline edges of cell-like structures, such as around 
(x,y) - (10", 40") (indicated by a dashed line). 



Each IBP is associated with a single patch, but a patch often 
contains multiple IBPs. 

The left-hand part of Fig. |4] shows the result of the patch 
analysis for the Can H sequence cutout shown in Fig. |2l The 
three panels display the corresponding x-y cutout, as well as 
the associated x-t and y-t slices through the binary cube made 
up of the 24 Ca ii H IBPs that are members of this particular 
patch. Each panel is integrated over x, y, or t, as appropriate, 
so that a larger blackness signifies the presence of more IBP 
pixels in the cube along the third dimension. Some IBPs visible 
in Fig. |2] are missing, such as the one around x x 4 arcsec 
and f 25 min, because they are not a member of this patch. 
The y-integrated x-t slice in the bottom-left panel shows the 
intersection at f w 12 min and the merger at f » 45 min that 
were already visible in the rightmost panels of Fig.|2] Adding 
the information in the y-t slice makes it clear that there are in 
fact three trails, with mergers at f « 20 min and f ^ 45 min. The 
apparent intersection in the x-t slice at f w 12 min is actually 
a disjunct crossing, but the merger at f ^ 45 min is real. The 
latter connects the two legs of the patch through the friends-of- 
friends labeling. 

Figure|4]also adds the corresponding triplet for the G band 
map sequence (three right-hand panels). There is a strong cor- 
relation between the two diagnostics, but Ca ii H indeed pro- 
vides more IBPs. The right-hand trail is very similar in the two 
y-t slices, but the left-hand trail mostly vanishes in the G-band 
sHce. 
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Fig. 6. Histograms of IBP lifetimes. Solid: G-band IBP life- 
times. Dotted: Ca ii H IBP lifetimes. IBPs with lifetimes shorter 
than 4 time steps (80 s) are rejected by our processing. 



3. Results and Discussion 

3.1. Patch patterns 

The resulting IBP collection consists of 387 features in the G- 
band map sequence that are classified as IPBs, and 848 IBPs in 
the Ca II H map sequence. There are more in the latter because 
IBPs are more easily identified in Can H. We find 149 G-band 
patches, 76 of which contain multiple IBPs, and 217 CanH 
patches, 125 of which contain multiple IBPs. A comparison 
shows that 85% of the G-band patches conincide spatially with 
Can H patches. 

The locations of all Ca ii H IBPs, irrespective of their time 
of appearance, are shown in Fig.|5] They show a striking pattern 
in which groups of IBPs appear to partially outline the edges of 
cell-like structures of mesogranular scale, e.g., around (x,y) = 
(10", 40"), where a large, conspicuous cell is marked with a 
dashed curve, around (20", 15"), and around (27", 42"). 

Similar patterns have be en reported before in studies of 



internetwork magnetog r ams (iDommguez Cerdena et al.ll2003 
Sanchez Almeida 2003 ; Dominguez Cerdeflal 2003h . In partic 



ular, iDomuiguez Cerdeflal t003 analyzed a 17-mm sequence 
of Fe I magnetograms and found that the stronger internetwork 
magnetic elements show a pattern that coincides with meso- 
granular upwellings. One may expect such a pattern to be set 
by underlying granular motions. Roudier & M uller (2004) ad- 
vected corks by measured granular flows and found that the 
corks concentrate at the boundaries of "trees of fragmenting 
gr anules". These were p reviously connected to mesogranules 
by iRoudier et al.l (l2003h . who found that such trees may have 
lifetimes of many hours. 

3.2. IBP lifetimes and number density 

Figure |6]shows a histogram of the durations over which our al- 
gorithm tracks the IBPs that are completely within the temporal 
and spatial boundaries, and that do not split or merge. We find 
an average IBP lifetime of 3.5 min in the G band and 4.3 min 
in Ca n H. Both the G-band and Ca ii H IBP lifetime distribu- 



tions show a long tail towards long Ufetimes, with maxima of 
19.3 min and 25.3 min, respectively. 

Berger et al.1 (Il998bh find longer lifetimes of 9.3 min for 



G-band network bright points. Possibly, their data and analysis 
permit better tracking of NBPs. Their method of NBP detection 
in G-band images by subtraction of a suitably scaled continuum 
image yields a much improved contrast between granulation 
and NBPs. However, there is no equivalent method for Can H 
images. In addition, internetwork bright points generally have 
substantially lower contrast than network bright points, and are 
therefore much harder to identify and track. By reducing the 
threshold value in our data reduction, IBPs mightbe followed 
through periods of lower intensity, which somewhat increases 
the mean lifetime at the price of many more misidentifications. 

Bright point lifetimes in internetwork regions may be 
shorter than in network regions. Where it is strong, the net- 
work disturbs the convection by fragmenting granules into "ab- 
normal granulation", so that flux tubes in network areas have 
a relatively quiet life compared to those in internetwork ar- 
eas, where full-size granules continuously crash into them. 
Their rapid, erratic movements make internetwork IBPs harder 
to track, and are furthermore likely to disturb the processes 
that make them bright. Since IBPs cluster in patches, we con- 
clude that, even though the IBP may momentarily become in- 
visible, the magnetic field element remains and may become 
bright agai n at some later t ime. T his agrees well with the con- 
clusion of iBerger & Title ( 200 lb that magnetism is a neces- 
sary but not a sufficient condition for the formation of a net- 
work bright point. We find an average IBP number density of 
0.02 Mm"^ in the G-band images and 0.05 Mm"^ in the Ca n H 
images. This is an order of magnitude less than the result of 
Sanchez Almeida et al. (2004), who find a number density of 
0.3 Mm"^ in their best G-band image. They manually located 
bright points in data with higher resolution than used in this 
analysis, which possibly allowed them to identify fainter fea- 
tures than our algorithm. Also, their 23" x 35" image contains 
a small network patch. We exclude all network and a fair area 
around it, w here the bright point number density is increased 
(cf. Fig. 2 of lSanchez Akneida et alll2004 . 



3.3. Statistical patch lifetime estimation 

Visual inspection of the patches identified by our algorithm 
shows that only a few patches begin or end within the dura- 
tion of the sequence, indicating that their typical lifetime sub- 
stantially exceeds one hour We therefore apply statistical argu- 
ments to obtain an estimate for the lifetime. 

Assume that we observe a patch with lifetime r in a se- 
quence lasting fobs- The patch is visible if the sequence starts 
between fobs before the patch emerges and the time that the 
patch disappears, i.e., in order to see the patch, the sequence 
must start in a time window of t -i- fobs- We distinguish three 
cases. 

(i) A patch persisting throughout the entire sequence. If the 
patch has a lifetime t > fobs, it is only visible for the entire 
duration of the sequence if it emerges up to t - fobs before the 
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start of the sequence. Therefore, the probability of finding such 
a patch is given by 



r 

Pi(t) — < T — fobs 
[ T + fobs 



if r < fobs , 

if T > fobs ■ 



(2) 



(ii) A patch that emerges as well as disappears within 
the duration of the sequence. A patch with lifetime t < fobs 
emerges as well as disappears within the duration of the se- 
quence if it emerges up to fobs - t after the sequence starts. The 
probability thus is 



^obs ^ 



if r < fobs : 



P2{t) - < T + fobs 

if T > fobs . 



(3) 



(iii) A patch that either emerges or disappears within the 
duration of the sequence. In case r > fobs, a patch is seen to 
emerge, but not to disappear, if the observation is started up to 
fobs before the patch emerges. Similarly, it is seen to disappear, 
but not to appear, if the observation is started up to fobs before 
the patch disappears. In the case that t < fobs, the patch is seen 
to emerge, but not to disappear, if the patch emerges up to r 
before the end of the observation, and is seen to disappear, but 
not to emerge, if the patch disappears up to t after the start 
of the observation. The probability pj, of seeing a patch with 
lifetime t either emerge or disappear in a sequence of duration 
fobs thus is 



P3(t) = 



2t 



T + fobs 

2 fobs 
T + fobs 



if T < fobs , 



if r > fobs . 



(4) 



We obviously have P\(t) + p2(j) + Pi(j) - 1. 

The expected numbers of patches A^' with / = 1,2,3 can 
be computed if we assume a lifetime distribution. We adopt a 
reasonable choice of an exponential distribution. 



-At 



(5) 



This distribution has one adjustable parameter. A, that is a mea- 
sure of the decay time scale of the distribution. The expected 
numbers of patches then follow from integrals of the form 



Jo 



N;=N PiiT)D,{T) dT, 



(6) 



where is the total number of patches observed. 

To obtain a fit for the parameter A, we need to count the 
number of patches of each type in our sequence. In nearly all 
cases a patch can be traced much longer by eye than that it 
is identified by the algorithm described in Sect. |2] We there- 
fore visually inspected the original data and identified those 
patches that remain visible during the whole sequence (A^i = 
124), those that both emerge and disappear during the se- 
quence (A^2 - and those that either emerge or disappear 
(A^3 = 68). We discarded the 14 patches that enter or exit the 
field of view during the sequence. We find an excellent fit for 
/I = 1.9 X 10"^ mm-\ yielding N[ = 123.7, N'2 = 8.6, and 



N'^ - 70.8. The excellent fit provides confidence in the valid- 
ity of the distribution D iij). The average patch lifetime is then 
given by 



<t) = I rD,,(T)dT= - s!530 + 50min, 
io ^ 



(7) 



where the error was estimated by variation of A^, within their 
error bounds set by Poisson statistics. 

4. Conclusion 

We have found numerous IBPs in photospheric and chromo- 
spheric quiet-sun internetwork cells. Our results show that 
Ca 11 H line-core filtergrams are well-suited for finding mag- 
netic flux tubes in the quiet-sun internetwork. IBPs in our 
Ca n H images show a good correlatio n with G-band IBPs. We 
follow ISanchez Almeida et al.l (120041) in attributing both IBPs 
in the G-band and Ca n H sequences to kiloGauss magnetic flux 
tubes. 

The IB P density that we measure is sig nificantly lower than 
the value of lSanchez Almeida et al.l (120041) . which is likely due 
to our lower resolution and more conservative identification 
method. 

We find that magnetic IBPs cluster into patches that are not 
homogeneously distributed over the internetwork, but rather 
seem to outline cell-like s tructures similar to the m a gnetic- 
element voids found by iDominguez Cerdena et alj (1200 3h 



marking mesogranular upwellings. This apparent mesogranu- 
lar distribution indicates that the magnetic elements which in- 
termittently appear as IBPs in a patch have a suflicient long 
lifetime to assemble at mesogranular vertexes. The strongest 
may eventually make it as NBPs to the supergranular bound- 
aries. Indeed, all IBP patches identified by our algorithm exist 
on time scales much larger than granular time scales. This re- 
sult seems to exclu de their a ttribution to a granular dynamo as 
described by, e.g., Cattaned ( 19991) . 

Through statistical analysis, we estimate an average life- 
time of about nine hours. This estimate depends on a distri- 
bution that admittedly cannot be verified from these observa- 
tions. This would require observations over much longer du- 
ration (> 10 hrs). The Solar-B mission in its sun-synchronous 
polar orbit will provide the high-resolution seeing-free long- 
duration observations that these analyses require. 
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